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Aims. For the first time we investigate the role of the grain surface chemistry in the Horsehead Photo-dissociation region (PDR). 
Methods. We performed deep observations of several H2CO rotational lines toward the PDR and its associated dense-core in the 
Horsehead nebula, where the dust is cold (Jdust - 20-30 K). We complemented these observations with a map of the p-H,CO 3o3 -2o2 
line at 218.2 GHz (with 12" angular resolution). We determine the H2CO abundances using a detailed radiative transfer analysis and 
. compare these results with PDR models that include either pure gas-phase chemistry or both gas-phase and grain surface chemistry. 

' Results. The H2CO abundances 2 - 3 x 10 '") with respect to H-nuclei are similar in the PDR and dense-core. In the dense-core 

the pure gas-phase chemistry model reproduces the observed H2CO abundance. Thus, surface processes do not contribute signifi- 
cantly to the gas-phase H2CO abundance in the core. In contrast, the formation of H2CO on the surface of dust grains and subsequent 
?^ , photo-desorption into the gas-phase are needed in the PDR to explain the observed gas-phase H2CO abundance, because the gas-phase 

chemistry alone does not produce enough H2CO. The assignments of different formation routes are strengthen by the different mea- 
^ ' sured ortho-to-para ratio of H2CO: the dense-core displays the equilibrium value (~ 3) while the PDR displays an out-of-equilibrium 

value (~ 2). 

Conclusions. Photo-desorption of H2CO ices is an efficient mechanism to release a significant amount of gas-phase H2CO into the 
C\| ■ Horsehead PDR. 

> ■ 

_^ ' Key words. Astrochemistry - ISM: clouds - ISM: molecules - ISM: individual objects: Horsehead nebula - Radiative transfer - 

' Radio Unes: ISM 

o 

^ ; 1. Introduction viouslv thought dOberg et alj|2009bllal: iMuiioz Caro et al.ll20Toh . 

1*—^ . These results led various groups to include photo-des orption into 

^ ■ Photo-dissociation region (PDR) models are used to understand PDR mo dels (see the results on H2O and O2 bv .Hollenbach et al.1 

the evolution of the far-UV illuminated matter both in our 2009: W alsh et al.ll20Tol; iHassel et al.ll20Toh . The availabihtv of 

Galaxy and in external galaxies. The spectacular instrumental vvell-defined observations is essential here to distinguish be- 

^> improvements, which happen in radioastronomy with the advent tween chemical assumptions about the significant grain sur- 

b<! °^ Herschel, ALMA and NOEMA, call for matching progresses face processes, i.e., adsorption, desorption, and diff^usion. It is 

in PDR modeling. In particulai-, the physics and chemisti-y of the now confirmed that some interstellar species are mosdy formed 

C3 dust grains and of the gas-phase are intricately intertwined. It is in the gas-phase (C O for instance), others on grains (CH3OH, 

■ ■ ■ ' well known that the formation of ice grain mantles leads to the Garrod e t al.1 12007|) . while the chemical routes for other com- 

removal of chemical reservoirs like CO, O, and other abundant piex species such as formaldehyde, are still debated because it is 

species from the gas phase, enabling new chemical routes to be ukely that soUd and gas-phase processes are both needed, 
opened and others to be closed. Despite their low temperature, 

the mantles are chemically active. Hydrogenation/deuteration re- Formaldehyde (H2CO) was the first organic molecu le dis- 
actions are known to be efficient, because hydrogen (or deu- covered in the interstellar medium (ISnvder et alJ Il969l) . It is 
terium atoms) can migrate on the surfaces, but reactions with a relatively simple organic molecule that can be formed in 
O, N, and C must also be considered. Complex molecules may the gas-phase and on the surface of dust grains. In the warm 
therefore be formed before they are released into the gas phase, gas, H2CO can trigger the formation of more complex organic 
Moreover, the release of the products into the gas phase hap- molecules dCharnlev et alj [1992). It is one of the most pop- 
pens either through thermal processes (evaporation) or non- ular molecules used for studying the physical conditions of 
thermal ones (cosmic ray or far-UV photon-induced desorp- the gas in astrophysical sources. Indeed, it js a good probe of 
tion). Recent photo-desorption experiments on water and CO the temperature and density of the gas ( iMangum & WoottenI 
ices show that this mechanism is much more efficient than pre- 1993i) . Owing to its large dipole moment (2.3 Debye), its 
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Table 1: Observation parameters for the maps shown in Fig.[T] The projection center of all maps is aiooo = 05''40'"54.27', 62000 = 
-02°28'00". 



Molecule 


Transition 


Frequency 


Instrument 


Beam 


PA 


Vel. Resol. 


Int. Time 


^sys 


Noise 


Obs.date 






GHz 




arcsec 




km s"' 


hours 


K(r;) 


K (r„b) 






Continuum at 1 .2 mm 




30m/MAMBO 


11.7X 11.7 















DCO+ 


3-2 


216.112582 


30m/HERA 


11.4X 11.4 





0.11 


1.5/2.0° 


230 


0.10 


2006 Mar 


p-HjCO 


3o3 - 2o2 


218.222190 


30m/HERA 


11.9X 11.9 





0.05 


2.1/3.4° 


280 


0.32 


2008 Jan 


HCO 


lo.i 3/2,2 -Oo.o 1/2,1 


86.670760 


PdBI/C&D 


6.69 X 4.39 


16 


0.20 


6.5* 


150 


0.09'' 


2006-2007 



° Two values are given for tlie integration time: tlie on-source time and tlie telescope time. * On-source time computed as if tlie source were always 
observed witli six antennae. ^ The noise values quoted here are the noises at the mosaic phase center (mosaic noise is inhomogeneous because of 
the primary beam correction; it steeply increases at the mosaic edges). 



Table 2: Observation parameters of the deep integrations of the 0-H2CO and P-H2CO lines toward the PDR and the dense-core. 



Position 


Molecule 


Transition 


y 


Line area 


Instrument 




6elT 


Int. Time 




RMS 


S/N 










[GHz] 


Kkms"' 








hours 


K (r„,b) 


K 






0-H2CO 


2l2 


-111 


140.839 


1.75+0.02 


30-m/C150 


0.93 


0.70 


1.9 


1.87 


0.061 


31 




p-HjCO 


2o2 


-loi 


145.603 


1.32+0.02 


30-m/D150 


0.93 


0.69 


1.9 


1.61 


0.047 


34 




0-H2CO 


2ii 


-lio 


150.498 


141 ±0.02 


30-m/C150 


0.93 


0.68 


14 


1.52 


0.053 


29 


PDR 


0-H2CO 


3l3 


-2,2 


211.211 


1.24+0.03 


30-m/B230 


0.91 


0.57 


1.1 


1.46 


0.096 


15 




P-H2CO 


3o3 


- 2o2 


218.222 


0.77+0.01 


30-m/B230 


0.91 


0.55 


3.9 


1.11 


0.052 


21 




p-HjCO 


3'>2 


-22, 


218.476 


0.17+0.01 


30-m/B230 


0.91 


0.55 


2.0 


0.27 


0.055 


5 




0-H2CO 


3l2 


-2„ 


225.698 


0.84+0.02 


30-m/A230 


0.91 


0.54 


6.5 


1.12 


0.079 


14 




0-H2CO 


2l2 


-1„ 


140.839 


2.56+0.01 


30-m/C150 


0.93 


0.70 


3.7 


3.46 


0.036 


96 




P-H2CO 


2o2 


-lo. 


145.603 


1.75+0.02 


30-m/D150 


0.93 


0.69 


1.9 


2.62 


0.044 


60 




0-H2CO 


2i, 


-1,0 


150.498 


1.89+0.01 


30-m/C150 


0.93 


0.68 


1.5 


2.52 


0.052 


49 




0-H2CO 


3l3 


-2,2 


211.211 


1.93+0.02 


30-m/B230 


0.91 


0.57 


2.0 


3.02 


0.065 


47 




P-H2CO 


3o3 


- 2o2 


218.222 


1.03+0.01 


30-m/B230 


0.91 


0.55 


3.0 


1.83 


0.057 


32 




P-H2CO 


3^2 


-22, 


218.476 


0.04+0.01 


30-m/B230 


0.91 


0.55 


4.5 


0.06 


0.037 


2 


Dense-core 


0-H2CO 


3,2 


-2„ 


225.698 


1.27+0.02 


30-m/A230 


0.91 


0.54 


84 


1.96 


0.073 


27 




o-H^^CO 


2l2 


-1„ 


137.450 


0.09+0.02 


30-m/D150 


0.95 


0.70 


2.0 


0.11 


0.063 


2 




p-HfCO 


2o2 


-lo. 


141.984 


0.10+0.01 


30-m/D150 


0.95 


0.70 


1.5 


0.11 


0.060 


2 




HDCO 


2ii 


-1,0 


134.285 


0.13+0.01 


30-m/C150 


0.94 


0.71 


2.0 


0.32 


0.042 


8 




HDCO 


3l2 


-2„ 


201.341 


0.05+0.01 


30-m/A230 


0.91 


0.59 


3.5 


0.13 


0.032 


4 




P-D2CO 


2l2 


-1„ 


110.838 


0.04+0.01 


30-m/AlOO 


0.95 


0.75 


4.9 


0.08 


0.031 


3 




0-D2CO 


4o4 


- 3o3 


231410 


0.04+0.01 


30-m/A230 


0.91 


0.53 


4.5 


0.09 


0.068 


1 



rotational lines are easy to detect from ground-based obser- 
vations. It is present in a variety of environm ents, such as 

proto-stellar 
young stellar 



Galactic HII regions ( e.g., [Downes et al 

cores (e.g., |Young et al 1l_2004l: iMaret et al 

objects (e.g., Araya et al. 2007), PDRs (e.g.,'Le urini et al.l20l6l) , 
starburst galaxies (e.g., Mangum et al..,.2008.) and comets (e.g., 
LSnvder etal.ll 19891: iJVIilam et al.ll2006h . 

The Horsehead PDR is particularly well-suited to investi- 
gate grain surface chemistr y in a UV irradiate d environment. It 
is viewed nearly edge-on (iHabart et al.l l2005h at a distance of 
400 pc (implying that 10" correspond to 0.02 pc). Thus, it is 
possible to study the interaction of far-UV radiation with dense 
interstellar clouds and the transition from warm to cold gas in 
a PDR with a simple geometry, very close to the prototypi- 
cal kind of source needed to serve as a reference to chemical 
models. Its relatively low UV illumination (j^^ = 60 in Draine 
units, [Draine 1978) and high density («h ~ 10^ - 10^ cm"^) 
implies low dust grain temperatures, from T^d ust ~ 30 K in the 
PDR t o Tdust ~ 20 K deeper inside the cloud ( Goicoechea et al.l 
l2009ah . The release of the grain mantle products into the gas 
phase is consequently dominated by photo-desoiption, while re- 
gions with warmer dust (the Orion bar or the star-forming cores) 
provide mix ed information on the thermal and non-thermal pro- 
cesses (e.g.. lBisschop et all2007h . 



In this paper we present deep observations of several 
formaldehyde lines toward two particular positions in the 
Horsehead neb ula: the PDR , corresponding to the peak of the 
HCO emission (iGerin et al. 12009), where the gas is warm (T\^ia ~ 
60 K); and the dense-core, a cold (Tkin < 20 K) condensation 
located less than 40" away from th e PDR edge, where HCO^ is 
highly deuterated (iPetv et al.l2007h . We present the observations 
and data reduction in Sect. |2] while the results and abundances 
are given in Sect. [3] In Sect. |4] we present the H2CO chemistry 
and PDR modeling. A discussion is given in Sect.|5]and we con- 
clude in Sect. |6] 



2. Observations and data reduction 

Tables [T] and |2] summarize the observation parameters for the 
data obtained with the IRAM-30m and PdBI telescopes. Fig. [T] 
displays the p-HjCO, HCO, DCO^ and 1.2 mm continuum 
maps. The P-H2CO line was mapped during 3.3 hours of good 
winter weather (~ 1 mm of water vapor) using the first polar- 
ization (i.e. nine of the eighteen available pixels) of the IRAM- 
30m/HERA single-sideband multi-beam receiver We used the 
frequency-switched, on-the-fly observing mode. We observed 
along and perpendicular to the direction of the exciting star in 
zigzags (i.e. + the lambda and beta scanning direction). The 
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H^CO 3q 3 2q 2 



[K.km/s] 



.Jy/beym] 




Fig. 1: Integrated intensity maps of the Horsehead edge. The 
intensities are expressed in the main-beam temperature scale. 
Maps were rotated by 14° counter-clockwise around the projec- 
tion center, located at (6x, Sy) - (20", 0"), to bring the exciting 
star direction in the horizontal direction and the horizontal zero 
was set at the PDR edge, delineated by the red vertical line. The 
crosses show the positions of the PDR (green) and the dense- 
core (blue), where deep integrations were performed at IRAM- 
30m (see Fig.O. The spatial resolution is plotted in the bottom 
left corner. Values of contour levels are shown on each image 
lookup table. The emission of all lines is integrated between 10. 1 
and 11.1 km s 



Table 3: Spectroscopic para meters of the obser ved lines obtained 
from the CDMS data base (iMiiller et al.ll200lh . 



IVTnlppiilp 


Transition 




E„ 


A , 




8ti 








[GHz] 


[K] 


[s-'l 






0-HtCO 


2,-. 


A 1 1 


140.839 


21.92 


5.3x10" 


-5 


15 


p-H,CO 




1 

loi 


145.603 


10.48 


7.8x10" 


-5 


5 


0-H2CO 


-^1 I 


- ill! 
^ lU 


150.498 


22.62 


6.5x10" 


-5 


15 


0-H2CO 


3 13 


-2i2 


211.211 


32.06 


2.3x10" 


-4 


21 


p-H,CO 


3o3 


- 2o2 


218.222 


20.96 


2.8x10" 


-4 


7 


p-H^CO 


322 


-22, 


218.476 


68.09 


1.6x10" 


-4 


7 


0-H2CO 


3l2 


-2n 


225.698 


33.45 


2.8x10" 


-4 


21 


o-H^'CO 


2l2 


-111 


137.450 


10.51 


4.9x10" 


-5 


15 


p-HfCO 


2o2 


-loi 


141.984 


2.37 


7.2x10" 


-5 


5 


HDCO 


2i, 


-lio 


134.285 


17.63 


4.6x10" 


-5 


5 


HDCO 


3l2 


-2ii 


201.341 


27.29 


2.0x10" 


-4 


7 


0-D2CO 


2l2 


-111 


110.838 


13.37 


2.6x10" 


-5 


5 


p-D,CO 


4o4 


— 3o3 


231.410 


27.88 


3.5x10" 


-4 


18 



p — H,CO 



10,5 X- 



21 

-^22 



O-H3CO 



33.5 
33.1 



-^12 
\3 



;2o 

■^21 



E, k (K) 



E.'k (K) 



22.6 
21.9 



15.4 
15.2 

E/k (K) 



multi-beam system was rotated by 9.6° with respect to the scan- 
ning direction. This ensured Nyquist sampling between the rows 
except at the edges of the map. The fully sampled part of the map 
covered a 150" x 150" portion of the sky. A detailed description 
of the HCO, DCO^ and 1.2 m m continuum obse rv ations and 
data re ductions can be found in 'Ger in et al] (|2009|) . iPetv et akl 
(l2007h . and|Hily-Blant et al. (2005). 

We performed deep integrations of 0-H2CO and P-H2CO 
low-energy rotational lines (see Fig. |2] and |3]l centered on the 
PDR and the dense-core. To obtain these deep integration spec- 
tra, we used the position-switching observing mode. The on-off 
cycle duration was 1 minute and the off-position offsets were (6 
RA, 6 Dec) - (-100",0"), i-e. the H 11 region ionized by crOri and 
free of molecular emission. From our knowledge of the IRAM- 
30m telescope we estimate the absolute position accuracy to be 
3". 

Th e data pro cessing was made with the GILDA^ soft- 
wares (|Petvll2005h . The IRAM-30m data wer e first calibrated to 
the Tl scale using the chopper-wheel method (Penz ias & BurrusI 
[Tom and finally converted to main-beam temperatures (Tmb) 
using the forward and main-beam efficiencies (Feff & Beif) 
displayed in Table |2] The resulting amplitude accuracy is ~ 
10%. Frequency-switched spectra were folded using the stan- 
dard shift-and-add method before baseline subtraction. The re- 



' See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the GILDAS softwares. 



Fig. 2: Lower energy rotational levels of para- (left) and ortho- 
H2CO (right). The energy above para ground-state is shown at 
the left of each level. The arrows indicate the transitions detected 
in the Horsehead. 



suiting spectra were finally gridded through convolution with a 
Gaussian to obtain the maps. 

3. Results 

3.1. H2CO spatial distribution 

The 218.2 GHz p-HjCO integrated line-intensity map is shown 
in Fig. [U together with the 86.7 GHz HCO, 216.1 GHz DCO+ 
integrated line-intensity maps and the 1.2 mm continuum- 
emission map. Formaldehyde emission is extended throughout 
the Horsehead with a relatively constant intensity. The H2CO 
spatial distribution ressembles the 1 .2 mm continuum emission: 
It follows the top of the famous Horsehead nebula from its front 
to its mane. It also delineates the throat of the Horsehead. The 
peak of the H2CO emission spatially coincides with the peak 
of the DCO^ emission, which arises from a cold dense-core. 
However, H2CO emission is also clearly present along the PDR, 
which is traced by the HCO emission. The PDR and dense-core, 
namely the peaks of the HCO and DCO^ emission are shown 
with green and blue crosses respectively. Gaussian fits of the 
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Fig. 3: Radiative transfer modeling of H2CO lines for two positions toward the Horsehead. Two left columns: the PDR position 
(rkin = 60 K, n(H2) = 6 X 10"^ cm"\ N(o-H2CO) = 7.2 x 10^^ cm"^) and two right columns: the dense-core position ( Tkin = 20 
K, n(H2) - 10^ cm"^, N(o-H2CO) = 9.6 x 10'^ cm"^). The two top rows display the ortho lines, for which we varied the column 
density around the best match (red curve) by a factor of 1.5 (blue curve) and 1/1.5 (green curve). The two bottom rows display the 
para lines, for which we kept the column density of the best match for 0-H2CO (red curves) constant and varied the ortho-to-para 
ratio of H2CO: o/p =1.5 (dashed blue), o/p=2 (dashed red) and o/p=3 (dashed green) 



H2CO lines at the HCO peak result in broader Une widths than at 
the DCO^ peak. That the lines are broader in the PDR confirms 
that H2CO lines toward the DCO^ peak arise from the dense- 
core rather than from the illuminated surface of the cloud. There 
is a peak in the H2CO emission toward the north-west region of 
the nebula, near the edge of the PDR, where two p rotostars have 
been identified (B33-1 and B33-28, Bowler et alj |2b09). These 
protostars heat the dust around them, so it is likely that H2CO 
has been evaporated from the grain ice mantles. 



3.2. H2CO column density 

We computed the column densities of H2CO at the PDR and the 
dense-core positions. For this we first used the Hj^CO lines to 
estimate the optical depth of the H2CO lines. Then, we made a 
first estimate of the column densities and excitation temperatures 
using rotational diagrams. Finally, we used these first estimates 
as an input for a detailed nonlocal non-LTE excitation and ra- 
diative transfer analysis to compute the H2CO abundances. The 



spectroscopic parameters for the detected transitions (shown in 
Fig. |2]i are given in Table [3] We assumed that the emission is 
extended and fills the 30m beam, as shown by the map of the 
3o3 - 2o2 transition (see Fig.lTJ. 

3.2.1 . Opacity of the H2CO lines 

We detected two transitions of the formaldehyde isotopologue 
Hj^CO in the dense-core position (see upper panels in Fig. |4|i. 
By comparing the flux between H2CO and Hj^CO for the same 
transition it is possible to estimate the opacity of the H2CO line, 
assuming that the Hj^CO line is optically thin, as follows: 

^ = (1) 

where /3 is the escape probability function, which in the case of 
a homogeneous slab of gas (de Jong et al. 1 980.) is equal to 

1 - exp(-3T) 



3t 



(2) 
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Fig. 4: Hj^CO and deuterated H2CO lines detected toward the 
dense-core. Gaussian fits are shown with red Unes. For HDCO 
and D2CO the Une width was fixed to the width of the HDCO 
(2n - 1 10) line, because it has the best signal-to-noise ratio. 



The isotopic abundance ratio '^C/'^'C - 60 (iLanger & PenziasI 
et al.ll2002l) is almost twice the line intensity ratio 
between formaldehyde and its isotopologue, and therefore the 
H2CO lines have moderate opacities. From the observations we 
estimate T2p-2,i ~ 1.6 and T2„,-i„, ~ 1.9 for H2CO in the dense- 
core. 

3.2.2. Rotational diagram analysis 

First-order estimates of the beam-averaged column densities and 
the rotational temperatures can be found by means of the widel y 
used rotational diagram analysis dGoldsmith & Langeii Il999l) . 
To do this, we assume that the gas is under LTE, and therefore 
all excitation temperatures are the same, and the energy levels 
are populated following Boltzmann's law. We built rotational di- 
agrams corrected for line-opacity effects through 



^thin ^ 

In-i^ -HlnC^ = In- 

J?« Z 



kTr„, 



(3) 



where is the total column density of the molecule, gu is the 
level degeneracy, Eu/k is the energy of the upper level in K, 
Z is the partition function at the rotational temperature T^ot, 
Ct = j^—;- < 1 is a line-opacity correction factor, where t is 
the opacity of the line, and A^*'" is the column density of the up- 
per level for an optically thin line when the source fills the beam. 
This last parameter is given by 



thin 



87Tkv^W 
hc^Aui ' 



(4) 



where k is the Boltzmann constant, v is the line frequency, W 
is the integrated line intensity, h is the Planck constant, c is the 
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Fig. 5: H2CO rotational diagrams corrected for line-opacity ef- 
fects at the PDR and dense-core position. Rotational tempera- 
tures are shown for each considered opacity. 



speed of light and A„/ is the Einstein coefficient for spontaneous 
emission. 

Ortho- and para forms of H2CO are treated as different 
species because radiative transitions between them are forbid- 
den. Resulting rotational diagrams are shown in Fig.|5]for three 
different 0-H2CO (2i2 - In) and p-HjCO (2o2 - loi) Hne- 
opacities (r = 0, 1 and 5). We find column densities of N ~ 
10'^ - 10'^ cm"^, depending on the opacity. We infer very dif- 
ferent rotational temperatures for 0-H2CO (Trot ~ 4 - 8 K) and 
P-H2CO (Tiot ~ 10-30 K), which are also lower than the well- 
known conditions in the PDR (rkin ~ 60 K) and in the dense- 
core (Tkin ~ 20 K). This suggests that the gas is far from ther- 
malization, and therefore we used these column densities and 
rotational temperatures as an input for a more complex analysis 
to derive the H2CO column densities. 

3.2.3. Radiative transfer models 

The critical density of a given collisional partner corresponds 
to the density at which the sum of spontaneous radiative de- 
excitation rates is equal to the sum of collisional de-excitation 
rates (7) of a given level 



ricriJK^K,, T'kin) - 



(5) 



Formaldehyde lines have high critical densities (~ 10^ cm"^, 
see Table HJ compared to the H2 density in the Horsehead 
(~ lO'* - 10^ cm"-'). Because we expect subthermal emission 
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Table 4: H2CO critical densities ( cm ) for three different col- 
liding partners computed for - 60 K. 





P-H2 


0-H2 


He 


2o2 


7.2 X 10" 


3.6 X 10' 


1.3 X 10" 


3o3 


1.6 X 10* 


9.9 X lO'' 


4.2 X 10* 


322 


5.8 X 10" 


4.7 X 10" 


2.5 X 10* 


2,2 


3.7 X 10" 


2.5 X 10" 


8.1 X 10" 


2,1 


4.3 X 10" 


2.2 X 10" 


8.7 X 10" 


3l3 


9.7 X lO'' 


7.0 X 10" 


2.3 X 10* 


3,2 


1.3 X 10* 


7.9 X 10^ 


3.2 X 10* 



{Tex ^ Tkin) for transitions with high critical densities com- 
pared to the H2 density, we used a nonlocal non-LTE radia- 
tive transfer code adapted to the H orsehead geometry to m odel 
the observed H2CO line intensities dGoicoechea et al.l2006 ). We 
used a nonlocal code to take into account the radiative coupling 
between different cloud positions that might affect the popula- 
tion of the energy levels. The code is able to predict the line pro- 
files. It takes into account line trapping, collisional excitation and 
radiative excitation by absorption of cosmic microwave back- 
ground and dust continuum photons. We included 40 rotational 
levels for 0-H2CO and 41 rotational levels for P-H2CO, where 
the maximum energy level lies at ~ 285 K for both species. We 
considered 0-H2, p-Hj and He as collision partners with the fol- 
lowing collisional excitation rates: 

- Collisional rates of 0-H2CO and p-HjCO with He are taken 
from Green (1991). 

- Collision al rates of 0- H7CO with 0-H2 and P-H2 from 
pTroscompt et alj ( l2009l) for the first 10 energy levels, i.e. 
Eu < 50 K. We compleme nted these data with He col- 
lision rates of iGreenI (Il991h scaled to H2. Following the 
new H2CO-H2 collisional rate calculations, we scaled the 
H2CO-He rates by a factor 2.5 instead of the usual ~ 1.4 
mass factor (A. Faure priv. communication). 

- Collisional rates of p-HjCO with 0-H2 and p-Hj from 
Troscompt et al., to be submitted. 

Results are presented in Fig. [3] for three different column 
densities. Best matches (see Table |5]l are for column densi- 
ties of N(o-H2CO) = 7.2 X 10^2 ^.^-2 N(p-H2CO) = 
3.6 X 10 '2 cm-2 in the PDR position, and N(o-H2CO) = 
9.6 X 10'2 cm-2 and N(p-H2CO) = 3.2 x lO'^ cm ^ in the dense- 
core position. In the excitation- and radiative transfer models we 
adopt an H2 ortho-to-para ratio of 3 (high-temperature limit), 
although it is li kely that the ortho- to-para ratio is lower in the 
Horsehead (e.g..'Habart et al.l l201lh . Indeed, the H2CO column 
densities are not sensitive to the change of the H2 ortho-to-para 
ratio for the physical conditions of the Horsehead (see Appendix 

3.3. H2CO ortho-to-para ratio 

The ratio between the column densities of 0-H2CO and P-H2CO 
provides information about the formation of the molecule, be- 
cause the characteristic conversion time from on e symmetry 
state to the other is longer than the H2CO lifetime dTudorie et all 
I2OO6 ). When the molecule forms in the gas-phase, a ratio of 3 
is expected, which corresponds to the statistical weight ratio be- 
tween the ground states of the ortho- and para species. A ratio 
lower than 3 is expected when the molecu le is formed on the 
surface of cold (rdus, $ 20 K) dust grains (iKahane et al.ll 19841: 
iDickens&Irvingi 19991) . From the derived column densities we 



Table 5: Column densities and abundances. 





Molecule 


rUK 


Dense-core 




Nil 


o . , 1 r\Z Z 

3.8 X 10 


6.4 X 10 




Nio-HiCO) 


7.2 X 10'2 


9.6 X 10'2 


Column densities 

[ cm--] 


Mp-HjCO) 
A'(HCO) " 


3.6 X 10'^ 
3.2 X 10" 


3.2 X 10'2 
<4.6x 10'^ 


iVI^rlUV^VJ ) 

yV(o-D2CO) 




1 A V 1 

5.1 X 10" 




Mp-DXO) I' 




5.1 X 10" 




[o-H,CO] 


1.9 X 10-'" 


1.5 X 10-'" 


Abundances 

rvi _ mx) 

' (JV(H)+2 W(H2)) 


[p-H,CO] 
[HCO] 
[HDCO] 
[D2CO] 


9.5 X 10-" 
8.4 X 10-'" 


5.0 X 10-" 
< 7.2 X 10-" 

2.5 X 10-" 

1.6 X 10-" 



" iGerin et aD( l2009h 

*> For Tex = 6 K (LTE). 



infer H2CO ortho-to-para ratios of ~ 2 in the PDR and of ~ 3 
in the dense-core. This suggests that in the dense-core H2CO is 
mainly formed in the gas-phase, whereas in the PDR H2CO is 
formed on the surface of dust grains. iDickens&Irvind (119991) 
measured H2CO ortho-to-para ratios between 1 .5 and 2 toward 
star-forming cor es with outflows , and r atios near 3 toward three 
quiescent cores. Hirgensen et al.l (12005 ) also found an ortho-to- 
para ratio of 1 .6 in the envelopes around low-mass protostars. 

3.4. HDCO and D2CO column densities 

We detected HDCO and D2CO in the dense-core (see two bot- 
tom rows in Fig.lUl, and we estimated their abundances assum- 
ing LTE. For Tex = 6 K we obtain N(HDCO) = 1.6 x lO'^ cm-^, 
N(D2CO) = 5.1 X 10" cm-2 and a D2CO ortho-to-para ratio of 
1, which translates into relative abundances or fractionation lev- 
els [HDCO]/[H2CO] = 0.11 and [D2CO]/[H2CO] = 0.04 for the 
inferred formaldehyde column densities in the dense-core. 

Deuterium fractionation can occur in the gas-phase by means 
of ion-molecule reactions, where D is transferred from HD 
to other species. High abundances of deuterated molecules 
compared to the elemental D/H abundance (~ 1.5 x 10-^, 
iLinskv et al.ll2006l) have been observed in different astrophys- 
ical environments, fr om cold dense cores and hot molecular 
cores even to PDRs. P etv et aP ('2007) found high deuteration 
([DCO+]/[HCO+] > 0.02) in the Horsehead dense-core. A pure 
gas-phase chemical model was able to repr oduce the observed 
fractionation level of HCO+ for Tkin < 20 K. lParise etan(l2009l) 
found high fractionation levels for DCN and HDCO toward the 
Orion Bar PDR ([XD]/[XH] ~ 0.01). They found that these ra- 
tios are consistent with pure gas-phase chemistry models where 
the gas is warm (> 50 K), so the deuterium chemistry is driven 
mainly by CH2D^, as opposed to colder regions (< 20 K) like the 
Horsehead dense-core, where H2D^ is the main actor Owing to 
the low temperature in the core it is likely that a non-negligible 
fraction of CO is frozen on the dust grains, enhancing the deu- 
terium fractionation. 

Another way to form deuterated molecules in cold environ- 
ments is trough D addition or H-D substitution reactions on the 
surface of dust grains (IHidakaetalJl2009l) . In the Horsehead 
core though, desorption from the grain mantles is not efficient 
in releasing products into the gas-phase (see sectionHJl. It is then 
more likely that the gas-phase HDCO and D2CO molecules de- 
tected here are formed in the gas-phase. Nevertheless, there can 
still be a considerable amount of deuterated H2CO trapped in the 
ices around dust grains. 
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Fig. 6: Photochemical model of the Horsehead PDR. Upper 
panel: PDR density profile (rtu = n(H) + 2n(H2) in cm"^). 
Middle panel: Predicted abundance (relative to n^) of H2CO 
(blue) and HCO (red). Lower panel: Predicted HCO/H2CO 
abundance ratio. In the two bottom panels, models shown as 
solid lines include pure gas-phase chemistry and models shown 
as dashed lines include gas-phase as well as grain surface chem- 
istry. The horizontal bars show the measured H2CO abundances 
and abundance ratios. 



4. H2CO chemistry 

We used a one-dimens i onal, steady-s tate photochemical model 
(iLe Bourlot et alJll993t lEe Petit et al.i.2006) to study the H2CO 
chemistry in the Horsehead. The physical conditions have al- 
ready been constrained by our previous observational studies 
and we keep the same assumptions for the density profile (dis- 
played in the upper panel of Fig. |6l), radiation field (x - 60) 
in D raine units), elernental g as-phase abundances (see Table 
6 in iGoicoechea et al.l l2009bl) and cosmic ray ionization rate 
(^ = 5 X 10-1^ s-i). 

Unlike other organic molecules like methanol, which 
can only be efficiently formed on t h e surface of grain s 
(iTielens & Whittetl [19971 Imoril I200I ICuppen et all l2009h . 
formaldehyde can be formed in both the gas-phase and on the 
surface of grains. Next, we investigate these two different sce- 
narios. 

4.1. Pure gas-phase chemistry models 

We use dthe Ohio State University (osu) pure gas-phase chemi- 
cal network upgraded to photochemical studies. We included the 
photo-dissociation of HCO and of H2CO (leading to CO and H2) 
with rates of 1.1x10'^ exp(-0.8Av) and 10"'* exp(-1.74Av) s"', 
respectively (Ivan Dishoe ck 1988). We also included the H2CO 
photo-dissociat i on channel that leads to HCO and H (see e.g., 
I Yin et al.ll2007t lTroell2007h with the same rate of the one that 



leads to CO and H2, and the atomic oxygen reaction with the 
methylene radical (CH2) to explain the high abundance of HCO 
in the PDR (Gerin et al. 2009). 

The predicted HCO and H2CO abundance profiles and the 
HCO/H2CO abundance ratio are shown as solid lines in Fig. |6] 
(middle and lower panel, respectively). The formation of H2CO 
in the PDR and dense-core is dominated by reactions between 
oxygen atoms and the methyl radical (CH3). The destruction of 
H2CO in the PDR is dominated by photo-dissociation, while it is 
dominated by reactions with ions in the dense-core. The pure-gas 
phase model satisfactorily reproduces the observed H2CO abun- 
dance in the dense-core {6x ~ 35 ) but it predicts an abundance 
in the PDR {6x ~ 15 ) that is at least one order of magnitude 
lower than the observed value. 

4.2. Grain chemistry models 

We considered the sur face chemistry reactions introduced by 
IStantcheva et al. I (l200l . which include the following sequence 
of hydrogen addition reactions on CO to form formaldehyde and 
methanol 

CO ^ HCO ^ H2CO HiCO CH3OH. 

H H H H 

We also introduce water formation via hydrogenation reactions 
ofO, OH until H2O. 

Adsorption, desorption and diff'usive reactions were intro- 
duced in the Meudon PDR code in the rate equations approach. 
The corresponding implementation will be described in a spe- 
cific paper (Le Bourlot et al., to be submitted) and we sim- 
ply mention the main processes included in the present study. 
We distinguish between mantle molecules, which may accu- 
mulate in several layers (e.g., H2O, H2CO, CH3OH), and light 
species (e.g., H, H2), which stay on the external layer Photo- 
desorption can be an eflicient mechanism to release molecules 
to the gas phase in regions exposed to s trong radiation field s, 
as shown recently in la boratory studies (lOberg et al.l l2009bllaL 
iMufioz Caro et al .1120101) . Thermal desorption is also introduced. 
It critically depends on the desorption barrier values, which are 
somewhat uncertain. Diffusive reactions occur on grain surfaces 
and the diffusion barriers are assumed to be 1/3 of the desorp- 
tion energy values. Photodesorption efficiencies have been mea- 
sured in the laboratory for CO, CO2, H2O and CH3OH. These 
experiments have shown that all common ices have photodes- 
orpti on yields of a few 10~^ m olecules per incident UV pho- 
ton (lOberg et al.l l2007i l2009allbllcb . Therefore, we also take a 
photo-desorption efficiency of 10"^ for those species that have 
not been studied in the laboratory. We assume in addition that 
for formaldehyde the two branching ratios toward H2CO and 
HCOh-H channels are identical, i.e. 5 x 10 Given the high 
density in the dense-core, the grains are assumed to be strongly 
coupled to the gas in the inner region, so that their temperatures 
become equal to 20 K in the dark region, whereas the illuminated 
dust grains reach temperature values of about 30 K. 

The predicted HCO and H2CO abundances are shown as 
dashed lines in Fig. |6] This model reproduces the observed 
H2CO abundance in the dense-core and predicts a similar abun- 
dance as the pure gas-phase model. This way, formation on grain 
surfaces does not contribute significantly to the observed gas- 
phase H2CO abundance in the dense-core. This is because of the 
low photo-desorption rates in the core caused by the shielding 
from the external UV field. On the other hand, the H2CO abun- 
dance can increase by up to three orders of magnitude in the 
illuminated part of the cloud (Ay < 4) when including the grain 
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surface reactions. The H2CO abundance now even peaks in the 
PDR, while it peaked in the dense-core in the pure gas-phase 
model. The model predicts a H2CO abundance peak in the PDR 
that is higher than the observed abundance averaged over the 
30m (~ 16")- This limited resolution prevents us from resolving 
the predicted abundance peak. Interferometric observations are 
needed to prove the existence of this peak in the PDR. 

5. Discussion 

H2CO has been detected in a variety of different astrophysical 
environments, with a wide range of gas temperatures and den- 
sities. It has been detected in diffuse clouds with high abun- 
dances (~ 10~^), observed in absorption again st bright HII re- 
gions (e.g.,' Liszt& Lucaslll995t [Liszt et al.l |2006). It is not well 
understood how H2CO can be formed and survive in such harsh 
environments, because gas-phase process cannot compete with 
the photo-dissociation and dust grain temp eratures are too high 
for molecules to freeze on their surfaces. iRoueff et al.1 (I2OO6I) 
detected absorption lines of H2CO at 3.6//m toward the high- 
mass protostar W33A, and estimated an H2CO abundance of 
~ I0~^ where t he gas has a temperature of ~ 100 K. Recently, 
iBergman etaP (l20Ilh found H2CO abundances ~ 5 x 10 in the 
p Ophiuchi A cloud core. Abundances of H2CO and other more 
complex molecules toward hot cores and protostars are high. In 
these regions the gas is dense and hot, so the dust grains also 
have high temperatures (> 100 K). Therefore, the ice mantles, 
formed in the cold pre-stellar phase, are completely evaporated. 
Once these molecules are in the gas-phase, they trigger an active 
chemistry in t he hot gas, forming even more complex molecules 
jCharnlev et al,.1992.) . 

H7 CO has also been observed in other PDRs. iLeurini et al.l 
(l20ld) detected H2CO in the Orion Bar PDR toward both the 
clump («H ~ 10'' cm"^) and the inter-clump («h ~ lO'* cm"^) gas 
components. They found higher H2CO abundances (~ 10 - 
10 than the ones inferred in this work for the Horsehead 
(~ 10"'°). Molecules trapped in the ice mantles can be thermally 
desorbed when the dust grains are warm enough. The dust tem- 
perature at which a significant amount of H2CO evaporates can 
be estimated by equating the flux of desorbing molecules from 
the ices to the flux of adsorbing molecules from the gas (see eq. 
5 in Hollenbach et al. 2009). Taking an H2CO desorption energy 
of 2050 K ( Garrod & Herb st 2006), we obtain an evaporation 
temperature of ~ 41 K. In the Orion Bar the dust grains have 
temperatures of Tdust > 55 - 70 K, so molecules can be desorbed 
from the icy mantles both thermally and non-thermally. But in 
the Horsehead PDR dust grains are colder (Tdust ~ 20 - 30 K), 
therefore molecules can only be desorbed non-thermally. Hence, 
the main desorption mechanism in the PDR is photo-desorption. 
In this respect, the Horsehead PDR offers a cleaner environment 
to isolate the role of FUV photo-desorption of ice mantles. In 
the Horsehead dense-core dust grains are also cold (~ 20 K), 
but photo-desorption is not efficient because the dust is shielded 
from the external UV field. Cosmic rays can desorb molecules 
from the ice mantles, but this contribution is not significant be- 
cause the desorption rates are too low compared to the H2CO 
formation rates in the gas-phase. Both the measured H2CO abun- 
dance and ortho-to-para ratio agree with the scenario in which 
H2CO in the dense-core is formed in the gas phase with no sig- 
nificant contribution from grain surface chemistry. 

We have shown that photo-desorption is an efficient mech- 
anism to form gas-phase H2CO in the Horsehead PDR. But, to 
understand the importance of grain surface chemistry over gas- 
phase chemistry in the formation of complex organic molecules. 



a similar analysis of other molecules, such as CH3OH and 
CH2CO, is needed. In particular, CH3OH is one of final products 
in the CO hydrogenation pathway on grain surfaces. It can also 
form H2CO when it is photo-dissociated. Therefore, their gas- 
phase abundance ratios will help us to constrain their dominant 
formation mechanism and the relative contributions of gas-phase 
and grain surface chemistry. Similar studies in different environ- 
ments will also bring additional information about the relative 
efficiencies of the different desorption mechanisms. 

6. Summary and conclusions 

We have presented deep observations of H2CO lines toward the 
Horsehead PDR and a shielded condensation less than 40" away 
from the PDR edge. We complemented these observations with 
a p-HjCO emission map. H2CO emission is extended through- 
out the Horsehead with a relatively constant intensity and resem- 
bles the 1 .2 mm dust continuum emission. H2CO beam-averaged 
abundances are similar 2 - 3 x 10 '*') in the PDR and dense- 
core positions. We infer an equilibrium H2CO ortho-to-para ra- 
tio of ~ 3 in the dense-core, while in the PDR we find a non- 
equilibrium value of ~ 2. 

For the first time we investigated the role of grain surface 
chemistry in our PDR models of the Horsehead. Pure gas-phase 
and grain surface chemistry models give similar results of the 
H2CO abundance in the dense-core, both consistent with the ob- 
servations. This way, the observed gas-phase H2CO in the core 
is formed mainly trough gas-phase reactions, with no significant 
contribution from surface process. In contrast, photo-desorption 
of H2CO ices from dust grains is needed to explain the observed 
H2CO gas-phase abundance in the PDR, because gas-phase 
chemistry alone does not produce enough H2CO. These different 
formation routes are consistent with the inferred H2CO ortho-to- 
para ratios. Thus, photo-desorption is an efficient mechanism to 
produce complex organic molecules in the PDR. Because the 
chemistries of H2CO and CH3OH are closely linked, we will 
continue this investigation in a next paper by studying the chem- 
istry of CH3OH in detail. 
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values for the H2 ortho-to-para ratio. We show models for an H2 
ortho-to-para ratio of 3 in red (high temperature limit), and we 
show models for the extreme case where the H2 ortho-to-para ra- 
tio is in green (low temperature limit). The difference between 
the models is less than 10%, which is within the observational 
uncertainties and therefore not significant. 
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Fig. A. 1 : Radiative-transfer modeling of H2CO lines for the core 
position in the Horsehead. The two top rows display the ortho 
lines and the bottom row displays the para lines. The best-match 
models are given in colors (Tkm - 20 K, n(H2) = 10^ cm"-', 
N(o-H2CO) = 9.6 X 10'2 cm-2, N(p-H2CO) = 3.2 x 10'^ cm^^), 
taking a H2 ortho-to-para ratio of 3 (red lines) and of (green 
lines). 



Appendix A: H2 ortho-to-para ratio 

We investigated the influence of the H2 ortho-to-para ratio 
adopted in the excitation and radiative transfer models. In 
Fig. lA. II we show the best-match models for the H2CO lines to- 
ward the core position in the Horsehead assuming two different 
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